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Abstract. We present mean-field calculations for the in-plane optical conductivity, the superfluid density, 
and the electronic Raman susceptibility in quasi two-dimensional systems possessing a ground state with 
two competing order parameters: d-wave density wave (dDW) and d- wave superconductor (dSC). In the 
coexisting dDW+dSC phase we calculate the frequency dependence of these correlation functions in the 
presence of impurity scattering in the unitary limit, relevant to zinc-doped cuprate superconductors. 
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The pseudogap phase of underdoped high-temperature su- 
perconductors (HTSC) has attracted considerable theo- 
retical 1 112131 115] and experimental }6I7I8I9] attention over 
recent years. Despite numerous efforts the microscopic ori- 
gin of the pseudogap and in particular its relationship with 
superconductivity still remains an open question. The de- 
struction of the Fermi surface above the superconducting 
transition temperature T c has been detected by a series of 
experimental techniques, both indirect and direct, includ- 
ing tunneling [91 and angle-resolved photoemission spec- 
troscopy [51715] . 

As to the theoretical description of this striking phe- 
nomenon, one class of models attributes the pseudogap 
phase as a precursor of superconductivity in which Cooper 
pairs form at a temperature T* but only acquire phase 
coherence at a lower temperature T c , where they form 
a uniform d-wave superconducting (dSC) condensate [2J. 
Within this framework the pseudogap and the supercon- 
ducting gap are intimately related. As opposed to this, 
another class of models invokes various phases which are 
not directly related to superconductivity, but rather com- 
pete with it. One of these is the d-wave density wave 
(dDW) scenario pQ, which has earned particular inter- 
est by now. The dDW is a prototype of an unconven- 
tional density wave (UDW), a particle-hole condensate 
with a fully developed gap of <i-wave symmetry. Its most 
interesting property is perhaps the fact that it lacks elec- 
tronic charge density modulation, hence the name "hid- 
den order" , and that there are real staggered orbital cur- 
rents circulating in the plaquettes of the Cu02 plane in 
an alternating fashion. Much attention has been devoted 
to show the consistency of the dDW scenario with sev- 
eral properties of HTSC. These include transport prop- 



E-mail: vanyolosOkapica . phy . bme . hu 



erties such as the in-plane and c-axis optical conductivi- 
ties [5110] , electronic Raman response [5] , thermodynamic 
properties [TT], and quasiparticle scattering on a local- 
ized impurity [12113114] . just to mention a few. Although 
the dDW scenario was introduced initially on phenomeno- 
logical grounds [T], its possible occurrence in microscopic 
models like the half-filled two-leg ladder has been checked 
since then [15] , 

As we have mentioned so far, the pure dDW phase has 
been studied quite extensively in recent years from many 
aspects due to its potential applicability to the anomalous 
normal state of underdoped cuprates. However, the pos- 
sible coexistence of this particle-hole condensate with d- 
wave superconductivity below T c and the consequences of 
the two-gap feature has not been fully developed, at least 
as far as the frequency dependent correlation functions are 
concerned [4116] , Mean- field calculations show that there 
is a substantial region on the temperature-doping phase 
diagram of the cuprates where the energetically stable so- 
lution of the gap equations is a mixture of dDW and dSC 
orders [4114117] . This coexisting dDW+dSC state has been 
studied within the context of the c-axis optical sum rule 
and superfluid density of HTSC [TU]. Also, there has been 
considerable interest devoted to the exploration of quasi- 
particle interference patterns and local density of states 
around a single impurity [12114115] . 

Partly motivated by recent infrared and electronic Ra- 
man scattering measurements showing evidence for two- 
gap feature in underdoped HTSC [19120] , the aim of this 
paper is to go beyond the study of single particle proper- 
ties of dDW+dSC and calculate two-particle Green's func- 
tions. Namely, we present analytical calculations of the 
in-plane optical conductivity and the electronic Raman 
susceptibility of dDW+dSC. Disorder effects are included 
within the unitary limit of impurity scattering relevant to 
Zn-doped high-T c cuprates [21122] . Also, coherence effects 
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are studied in detail: we find that while electromagnetic 
absorption (conductivity) is a clear type I process in pure 
(single component) dDW, electronic Raman response can 
be either type I or II, depending on the symmetry of the 
Raman vertex. In particular we find that in Ai g and B\ g 
channels the dDW coherence factor is type I, whereas in- 
terestingly in B 2g it changes to type II. This Raman exam- 
ple shows that one has to be rather careful in determining 
the coherence factor of a response function in the differ- 
ent condensates, i.e. in dDW and dSC, because the general 
rule of thumb according to which one only has to change 
the coherence factor from type I to II (or vice versa) when 
switching between dDW and dSC, might not work every 
time. It works for example in the case of infrared response, 
but not necessarily in the Raman experiment. Further ex- 
amples are given in Section [31 

The paper is organized as follows. In Sec. [5] we describe 
our model and the analytic formalism. In Sec. [3] we calcu- 
late the optical conductivity and the superfluid density in 
dDW+dSC, while in Sec. U] electronic Raman scattering is 
studied in the relevant A\ g , B\ gi and B 2g symmetries. Fi- 
nally, in Sec. [S]we conclude and summarize. An Appendix 
is also given for some technical details. 



2 Formalism 

Let us consider the following effective low-energy Hamil- 
tonian |23l24l25] 

H o = ^(e(k) ~ ^) c k,cr c k,<7 

k,cr 

+ ^ (^l(k)Ck )0 .Ck+Q,<r + h.c. 



k.cr 



+ E(^( k )<t c -k,4 + h - C 



(1) 



where e(k) is the electronic kinetic energy, \x is the chem- 
ical potential, Q = (w/a,Tt/a) is the dDW nesting vector, 
and the amplitudes and momentum dependent parts of 
the dDW and dSC order parameters, respectively, sepa- 
rate via A 1 (k) = A 1 f(k) and A 2 (k) = A 2 f(k), with 



/(k) = (cos(afca;) — cos(afc y ))/2. 



(2) 



Now the Green's function of the pure system without im- 
purities is obtained as [24125] 

G (k, iwn)" 1 = iui n - £(k)p 3 <7 3 + rj(k)a 3 

+ A 1 f(k)p 2 - \A 2 \f{k)e^p 3 a u (3) 

where £(k) = (e(k) — e(k + Q))/2, the imperfect nesting 
term [23124125] reads 



7 ? (k)= / i-(e(k) + e (k + Q))/2, 



(4) 



and the corresponding four-dimensional spinor field is de- 
fined as 

<F+(k) = (c+ t , c_ k ,+, c£ +Q)t , c_k-cu). ( 5 ) 



As is well known from dDW theory [115118] , the order pa- 
rameter of dDW is purely imaginary. Thus, by factoring 
out the imaginary unit A\ must be necessarily real. To 
keep the following calculations as simple as possible, but 
without the loss of generality we may choose the dSC or- 
der parameter A 2 to be real too, because its phase (f> 2 is 
unrestricted in our mean-field theory. From the poles of 
Go(k, uS) the quasiparticle excitation spectrum is obtained 
as uj = ±E±(k), with 



25±(k) = y(±^(k)+4?(k)-»,(k)J +Al(k). (6) 

Impurity scattering can be incorporated into the the- 
ory along the lines developed first by Abrikosov and Gor- 
kov in the early sixties [26127128] . The interaction of elec- 
trons with non-magnetic point-like impurities is described 
by the Hamiltonian [29130131] 



Hi = \ ^ lqR ^ + (k + qMR,Mk), 



(7) 



k,qj 



where the potential in the four-dimensional Nambu-Gorkov 
space reads 



U 



( 1 


V o 





-1 





1 

_ e -iQRj o 



\ 



-1 / 



(8) 



with Uj being the position of the jth impurity atom and 
U is the amplitude of the delta-potential. Following the 
method of Ref. [29] , the electronic self-energy from impu- 
rities is 



S R {iu> n ) = m (uiH)- 1 - ^G(k,^)j 



(9a) 



where rij = N\/V is the impurity concentration and the 
subscript R means the position of an impurity over which 
the average should be taken. Its diagrammatic represen- 
tation is shown in Fig. [TJ But most importantly, because 
of the requirement of self-consistency, G(k, iu> n ) is itself 
the renormalized Green's function appearing in Dyson's 
equation 



G(k, iuj. n 



G (k, iu> r , 



i 



(9b) 



Before proceeding further with the analytical calcu- 
lation we have to comment on an important simplifica- 
tion we will employ. The momentum dependence of the 
quasiparticle band structure given by Eq. ^ is very com- 
plicated because of the finite imperfect nesting term. As 
is well known by now from ARPES measurements un- 
derdoped cuprates do not exhibit perfect nesting and the 
experimental situation is best described by a t-t' model, 
where both the next-nearest neighbor hopping term t' and 
the doping fi is finite [16118] . This gives rise to finite rj, and 
consequently to a rich temperature-doping phase diagram 
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Fig. 1. Infinite series of non-crossing self-energy (or T- 
matrix) diagrams due to impurity scattering. The solid 
line denotes the electron propagator G, while the dashed 
line is for electron-impurity interaction. Dashed lines com- 
ing from the same cross represent successive scattering on 
the same impurity. 



[14117] . However, in our model calculation focusing on the 
coherence effects coming from the different components of 
the dDW+dSC condensate, in a first approximation it is 
sufficient to study the r\ — case only. As we shall see 
shortly, the effect of coherence factors on the lineshapes 
of the dynamical susceptibilities are rather robust and re- 
main strong in a system with imperfect nesting. An inclu- 
sion of a finite but small rj (k) would of course introduce 
additional features in the spectra, but this would require a 
fully numerical approach and is therefore out of the scope 
of the present investigation. Nevertheless, further devel- 
opments in this direction are highly desired. 

Neglecting imperfect nesting we find that electron or 
hole pockets will not form around the dDW gap nodes and 
the two gap functions add up squared as the quasiparticle 
energies become doubly degenerate 



£±(k) = E(k) = y/?M + (Aj + Al)P(k). (10) 

It is clear that these quasiparticles behave in some respect 
the same as those in dSC. Now, however, they feel a com- 
bined effective d-wave gap A(k) = Af(k), where 



A 



Al 



(11) 



where K(z) is the complete elliptic integral of the first 
kind, a = r/A is the dimensionless pair-breaking pa- 
rameter, r — 2ni/(-7riVo) is the scattering rate, and iVo 
is the normal state density of states per spin. In arriving 
to Eq. (1131) the momentum integrals were performed ana- 
lytically in the continuum limit, where the Fermi surface 
is assumed to be cylindrical, the d-wave angular depen- 
dence simplifies to /(k) = cos(2</>), and thus the Fermi 
surface average is simply an average over <j> [32 33 39 . We 
note that this commonly used simplification of momen- 
tum integrals is purely technical. It allows for analytical 
calculations, but at the same time captures the essence of 
physics in the sense that the qualitative features remain 
the same as those obtained from numerical integration 
with the exact lattice periodic functions |5I40] . The renor- 
malized Matsubara frequency at zero external frequency, 
namely Cq = u)(lo„ = 0), plays a central role in the theory 

gnu. it is given by 



This mapping between a dSC and the particle-hole sym- 
metric dDW+dSC involves that the whole thermodynamic 
behavior of the latter, as well as its response to static 
external perturbations are the same as those of a dSC 
|32l33j . because the density of states (DOS) is itself the 
same. Only the energy scale is different according to Eq. ([TTj) 
By direct substitution one can also verify that in unitary 
limit, where U — > oo, the solution of Eqs. © is 



G(k, w n ) 1 = iid.„ 



£(k)p 3 o- 3 

+ A 1 f(k)p 2 -A 2 f(k)p 3 a 1 , (12) 



where the Matsubara frequency is the only quantity to be 
renormalized by impurities. The detailed analysis of the 
unitary limit is covered in many textbooks 32 33 34 35 36 
and publications [29 30137138] , therefore we shall not re- 
peat it here. We shall only summarize some essential ex- 
pressions that will be utilized repeatedly during the cal- 
culation of the complex conductivity and the Raman re- 
sponse. The analytically continued renormalized frequency 
u = iuj(ioj n — > lo + iO)/A obeys 



(13) 



K 



1 



(14) 



Apropos of Eq. (ITT|) we have already pointed out that 
the single particle properties of a dDW+dSC state with 
rj set to zero are the same as those of a dSC due to the 
matching DOS |32I41] . The similarities, however, end here 
because at finite frequency the dynamical response is dom- 
inated by coherence effects. These reveal the true nature 
of the ground state, the fact that it is actually a competing 
superposition of dDW and dSC condensates, both with its 
own order parameter. 



3 Optical conductivity 

In this section we calculate the frequency dependence of 
the in-plane optical conductivity and the superconducting- 
condensate fraction at zero temperature for the coexist- 
ing dDW+dSC phase with impurities. Following standard 
methods [27135143] the complex conductivity is obtained 
from the correlation function x = ([Ja>Jp]), where the 
paramagnetic current operator reads in our four-dimen- 
sional Nambu notation as 

J a (q) = -e^ Q (k + q/2)<F+(k)p 3 ^(k + q). (15) 

k 

In this expression it is the p 3 matrix structure that should 
be really kept in mind, because this is the essential ingre- 
dient here that distinguishes the infrared response from 
the Raman response via the different coherence factors. 
As to the specific momentum dependence of the quasipar- 
ticle velocity, in accord with our assumption of a cylindri- 
cal Fermi surface, we use u a (k) = fc Q /m, with m being 
the effective mass of the particles. With this the current 
correlator is 



X{iv n ) = S a /3 



47re 2 n T 
m 2 V 

-^ + e 2 (k)-^(k) + ^i(k) 

(C»+i52(k))(fl«+25»(k)) 



(16) 
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where n is the total electron density, and the prime on uj' n 
means that the renormalizcd Matsubara frequency has to 
be computed at the shifted frequency u) n + v n . 

At this point it is worth stopping for a moment and 
discuss coherence effects in the different condensates. In 
BCS theory [36 44 the coherence factors show up in the 
calculation of two-particle Green's functions, a typical ex- 
ample of which is the current-current correlator we are 
just studying in Eq. (fT6"|). In pure superconductors the 
frequency sum can be performed 43 45 , leading to the 
appearance of the coherence factors. For brevity we only 
quote here the long wavelength result associated with the 
creation of two quasiparticles, it is: 1 — (£ 2 =p A 2 )/E 2 . The 
upper and lower signs refer to type I and type II coher- 
ence factors, respectively. Electromagnetic absorption is a 
representative of type II processes, where the factor van- 
ishes exactly, signaling the well known fact that at finite 
frequency the optical response of a pure superconductor is 
zeroQ Coming back to our original formula in Eq. (j!6[) we 
see that due to the presence of impurities the Matsubara 
sum cannot be done analytically anymore, instead, it is the 
momentum integration that will be performed first. Nev- 
ertheless, it is not too difficult to figure out even from this 
expression that the plus sign in front of the superconduct- 
ing order parameter would involve clear type II coherence 
effects if the ground state was a single component dSC. 
As opposed to this in a single component dDW, where 
only ^\i(k) exists, the sign is just the opposite, leading to 
clear type I behavior. From this example one might have 
the impression that the same physical quantity is always 
of different type in a superconductor and a density wave. 
Another example of great importance is the charge den- 
sity with (T3 matrix structure, which has plus and minus 
signs in front of the dDW and dSC order parameters, re- 
spectively, i.e. the situation is just the opposite compared 
to the charge current. This general rule of thumb seems to 
work usually, and is related to the different structures of 
the Nambu space in these condensates. In density waves 
it is the k > k + Q transformation that spans the space, 
and the spinor & consists of particle type operators only. 
In superconductors, on the other hand, the transformation 
properties of the physical quantities are studied under in- 
version in momentum space, that is under k — > — k, and 
the spinor mixes particle and hole type operators. This 
extra canonical transformation from particles to holes in- 
troduces an additional minus sign compared to that in 
a density wave. But there are interesting exceptions as 
well. One of these is the electronic Raman susceptibility, 
which, as we shall see in Section |4j turns out to be a type 
I process in superconductors. In dDW, however, the sym- 
metry of the Raman vertex decides whether it is type I or 
type II. Another exception is the spin density whose ma- 
trix structure is simply the identity, leading to plus signs 
in the numerator and type II behavior in both conden- 



1 Note, however, that though the real part of the conduc- 
tivity vanishes at every finite frequency, the conductivity sum 
rule requires the optical weight to be conserved, which means 
that all spectral weight is necessarily shifted to uj — in the 
form of a Drude peak (DS(uj)). 



Physical quantity Nambu structure dDW dSC 

Charge current pz — + 

Charge density 0-3 + 

Spin density 1 + + 

Raman A\ B P303 — — 

Raman B\ g pzoz 

Raman B29 03 + — 

Table 1. Coherence factors associated with a few physi- 
cal quantities in dDW and dSC condensates. The ± signs 
appear in front of the order parameters in the one bub- 
ble calculation of the relevant correlation functions. For 
examples see Eqs. p^|) and (|3"Tj) . The minus (plus) sign 
is commonly referred to as a type I (II) coherence factor 



sates. In order to give a brief and easy-to-use reference of 
all these results we gathered these physical quantities and 
the appropriate coherence factors in Table [1] with special 
emphasis on the difference between the dDW and dSC 
cases. 

Now, after carrying out the momentum integration in 
Eq. (|16p and transforming the frequency sum into a con- 
tour integral the (regular) real part of the optical conduc- 
tivity can be expressed as [29130141] 



(jAUj) 



2 n I n (uj) + I ph (uj) 



mix A uj 
where, using density wave terminology 



(17) 



In(w) 



da; ^tanh [^^- 



tanh ( — — ) 
\2TJ 



x Re |f(m(x + uj), u(x)) - F(u(x + uj), u(x))\ (18e 
is the so-called normal (intraband), while 

x Re {f(u{x - uj), u(x)) ~ F(u(x - u),u(x))\ (18b) 

is the pair-breaking (interband) contribution. Further, the 
dimensionless function F is 



dxtanh ( — ] 

^t) 



F(u,u') = 



2^1 
A 2 



1 



u' 2 — u 2 



— u(u + v!) 



K 



VI - u 2 E + 



u 2 K 
Vl -u 



M It 



(19) 



where the argument of K and E, the complete elliptic in- 
tegrals of the first and second kind, is 1/ yl — u 2 . This is 
the generalization of the corresponding F functions of dSC 
[H] and UDW [33]. Strictly speaking, as it stands above, 
it yields only the quasiparticle contribution to the conduc- 
tivity. However, due to the fact that in our approximation 
the scattering potential U (see Eq. ©) is momentum in- 
dependent, the vertex corrections exactly vanish and thus 
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Eq. ([TO)) gives the total contribution to <j\. This is not 
necessarily so in the Raman calculation, as we shall see 
shortly in Section 2] because there the ladder diagrams 
can renormalize the quasiparticle result, for example in 
the Big symmetry, and for this to see the vertex equation 
has to be solved too. 

The conductivity at zero temperature in the coexisting 
dDW+dSC phase is shown in Fig. [2] for different values of 
impurity concentrations. Note that on different curves the 
only parameter varying is a. The gap amplitudes are fixed, 
with the constraint in Eq. (jlip of course, and their renor- 
malization due to impurity scattering is not taken here 
into account^ Rather, we can think of these as experi- 
mentally adjustable parameters of the model. We can see 
that the stronger the pair-breaking is, the more the line- 
shapes approach the metallic Lorentzian form. One could 
think, naively, that the robust peak around U) = 2A must 
be the result of the density of states, in particular its Van 
Hove singularities located at the gap edge uj = ±A |33I41] . 
As the direct calculation verifies, this argumentation holds 
only in unconventional density waves [29130] , like the dDW 
itself [5], but in dSC the type II coherence factor sup- 
presses this transition [H]. Thus, the fact that this strong 
peak still develops, can be attributed to the dDW compo- 
nent of the ground state alone. The small bump at uj = A 
results from transitions from and to the narrow impurity 
band around the Fermi energy [?2 30 . It is also worth 
mentioning that in the dc conductivity the a — > and 
uj — > limits cannot be interchanged, this is why we ob- 
tain different result in the pure case. We believe that the 
right procedure is the one where the uj — > limit is taken 
first [30]. 

The dc conductivity in the ground state is obtained 
from Eq. ([17]) as 

0,(0) = G \ - 1 E ( -=L=) , (20) 

where Co has been defined previously in Eq. (1141) . This 
is formally the same as the one in dSC [H], but again 
A is the amplitude of the effective gap, which is always 
larger than that of the dSC component, A2. As the terms 
containing Co on the right hand side are always of order 
unity, and practically do not depend on impurities, this 
result for the dc conductivity is a manifestation of Lee's 
universality relation [47] . 

It might be of interest to examine the pure case ana- 
lytically too. For a = we find 

, . e 2 n f Ai\ 2 , / uj \ 

ai{w) = ^A\-A) tanh b0 

x t L (x- 1 (K(x)-E(x)) x<l, 

2xyK{x~ 1 )-E{x~ 1 ) x > 1, 1 ' 



2 For this we would have to solve the coupled gap equations 
for Ai and A2 self-consistently, for each a, but this enterprise 
is out of the scope of this publication. 




uj/A 



Fig. 2. (Color online) Zero temperature optical conduc- 
tivity of dDW+dSC in the unitary limit for a = 0, 0.01, 
0.025, 0.05, 0.1, 0.2, and 0.5. The curve belonging to a = 
meets the vertical axis at u = 3 A at the lowest. Further 
A t / A = 2/3. 

where x = uj/(2A). Perhaps this is the result that shows 
most concisely what we have discussed above apropos of 
Fig. [5] Namely, that it is the dDW component of the 
dDW+dSC condensate that leads ultimately to a finite 
optical response, and in particular to the logarithmic sin- 
gularity at the maximal optical gap 2 A. If the ground 
state consists entirely of a single component dDW, with 
A = Ai, this expression correctly simplifies to the UDW 
result |ig] . 

The fact that even in the pure case the dDW+dSC 
state exhibits finite absorption implies the weight of the 
Drude peak should decrease. In other words due to the 
presence of the dDW component spectral weight must be 
transferred from zero to finite frequency. This decrease of 
the Drude weight can be summarized nicely in the form 
of the in-plane superfluid density p s (T, r), which is calcu- 
lated at T = as 

The a rgumen t of the complete elliptic integrals K and E 
is 1/vl + u 2 . It is very interesting to observe that what 
appears in the large parentheses is formally the conden- 
sate fraction of dSC [42] . More remarkable is, however, 
that because of the prefactor not even in the pure sys- 
tem would the superfluid density reach unity [23124) . The 
missing electrons are in fact not missing but are frozen in 
the dDW component of the ground state. Recall that in 
a density wave, unless the phase fluctuations of the order 
parameter are taken into account, there is no dc conduc- 
tivity in the ground state [49150] . 
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4 Electronic Raman scattering 

In this section we proceed with the calculation of the elec- 
tronic Raman susceptibility in dDW+dSC with impuri- 
ties. The calculation is done along the same lines as the 
complex conductivity. In this case, however, instead of the 
current, the correlation function of the effective density 
has to be computed, where the velocity vertex is replaced 
by the Raman vertex 7(k) [51152) . The operator of the 
effective density is 

p(q) = ]>>( k + q/ 2 ) c iU c k+q,- (23) 

k,cr 

In the long wavelength limit the vertex is related to the 
generalized mass tensor of electrons [40153154] . We, how- 
ever, instead of making an ansatz for its momentum de- 
pendence, which obviously depends strongly on the details 
of the dispersion e(k), expand it as usual with the basis 
functions of the tetragonal point group D±h relevant to 
cuprate superconductors 39 55 56 

7 (k) = £ 7 £A£(k). (24) 

Here fi represents the even parity irreducible representa- 
tions of the group, L stands for the Lth-order contribu- 
tion to the vertex which transforms according to the /zth 
irreducible representation, and the unknown coefficients, 
responsible for the overall scale of absorption in the given 
channel, will be thought of as free parameters to be fit- 
ted to experimental data. As to the basis functions the 
following low order contributions will be considered [5I40J 

{l,cos(afc :E )-|-cos(afc y ) fi = Ai g , L = 0,4, 
cos(ak x ) — cos(aky) — Bi g , L = 2, 

sin(afc x ) s'm(aky) [i = B2 g , L = 2. 

(25) 

It is very interesting to see that among these examples 
neither the L = term in the A\ g symmetry, which is 
apparently nothing else than the density vertex, nor the 
Bi g basis function changes sign under the transformation 
k — > k + Q. This in turn leads to the observation that 
in Nambu space these contributions will exhibit a matrix 
structure equivalent to the charge density, that is 173. In 
stark contrast to this, all other basis functions do change 
sign leading to a distinct /O3173 behavior, and as we shall 
see in Eq. (|3"Tj) . to different coherence factor as far as the 
dDW is concerned. With these, our task is now reduced 
to finding the correlation functions \l = ([Pl^Pl\) °f the 
channel dependent operators 

p£(q) = iiY,K(w + mmk+<D, (26) 

k 

where the Nambu matrix M is either 173 or P3CT3, depend- 
ing on the case considered. 

Let us start the derivation with the definition of the 
correlation function we shall compute |52) . At long wave- 




Fig. 3. Diagrammatic representation of the vertex equa- 
tion in Eq. ([2"5f . The dot is the bare vertex, the shaded 
triangle represents the vertex correction due to impuri- 
ties, the cross plays the role of an impurity, the solid line 
is the Green's function, and the double dashed line denotes 
the self-consistent T-matrix, whose diagrams are shown in 

Fig.rrj 

length and imaginary frequency we find 

x l(iv n ) = ~( 7 D 2 | E A£(k)Tr(A/G(k^„) 

x Al(k,iuj n ,iuj n + iis n )G(k,iuj n + iv n )), (27) 

where G is the Green's function from Eq. (TT2l , and for our 
approximation to be conserving [43] . the impurity renor- 
malized vertex A^ L has to obey the following vertex equa- 
tion 

yl£(k, ioj n , iuj n + iv n ) = A£(k)M 

+ y^E j^^ T R( iuJ n)G(p,iuj n )Al(p,iuj n ,iui n +iv n ) 

p R 

x G(p, iuj n + iis n )T R (iuj n + iv n ). (28) 

Its diagrammatic representation is shown in Fig. [3J In 
the full non-crossing approximation the T-matrix incor- 
porates the effect of impurity potential to infinite order. 
Its diagrams are the same as those of the self-energy in 
Fig. [TJ Correspondingly, its analytical expression is given 
by Eq. (|9"aj) too, only the factor ri\ needs to be ignored. In 
the unitary limit, where U — > 00, it simplifies to 

Now, making use of these results and the fluctuation- 
dissipation theorem, the Raman intensity is found to be 
proportional to the imaginary part of the retarded corre- 
lator 

Imx£M = - 2N ^} Y (W + JphH) , (30) 

where the normal and pair-breaking integrals are given 
by Eqs. (fT8|) . but now the F function becomes channel 
dependent and needs to be calculated for every symmetry 
separately^ These cases will be studied in the following 
subsections. 

3 The conductivity and Raman calculations suggest that 
within mean-field theory every correlation function can be ex- 
pressed with these two integrals in the long wavelength limit. 
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Before turning our attention to the channel dependent 
spectra and the full solution of the vertex equation, it 
is worth writing down the quasiparticle contribution to 
the susceptibility. Diagrammatically this is the one bub- 
ble result, the only term that survives in a pure system, 
and its analytical expression is very similar to that of the 
current-current correlator in Eq. (|16p . It is obtained by 
approximating the full vertex by the bare one 



Xi(*fn)qp = - 

' ^A£(k) 2 



4(7D 2 



-u ln o J ' n +e(k)±Ai(k) 



(w*+^(k))(«« + ^(k)) 



(31) M 



where the plus and minus signs in front of the dDW gap 
belong to M — cr 3 and M — /O3O3, respectively. From 
this we can see explicitly what we have already indicated 
before, that (i) the appearance of 03 leads indeed to a sign 
change in front of the superconducting gap compared to 
Eq. (|T6|) . turning the Raman scattering to a type I process 
in superconductors, and (ii) in dDW it depends on the 
symmetry of the vertex whether inelastic light scattering 
is a type I or II process, see Table [TJ 



4.1 Raman spectra in A\ g symmetry 

In the fully symmetric representation the leading terms in 
the momentum dependence of the Raman vertex are given 
by Eq. (j2"5)l as 70 + 74A4(k), where for brevity we drop the 



superscript fi = A\ g . As pointed out before, the first term 
is associated with 03 in Nambu space, which means that 
its contribution to the full A\ g response is essentially the 
correlator of the charge density. However, due to the long 
range Coulomb interaction between electrons light scatter- 
ing on intercell charge fluctuations in the long wavelength 
limit is totally screened out, thus its contribution vanishes 
exactly [2D] . 

The L = 4 term does not couple to long wavelength 
density fluctuations, instead it is responsible for intracell 
charge fluctuations. Therefore, the total contribution to 
A\ g Raman intensity comes from this term. First we have 
to address the issue of vertex corrections. For this we 
have to solve Eq. (j28l) with M — ,03 03 and Az,=4(k) = 
cos(afc;;) + cos(afcj,). The solution turns out to be rather 
simple: vertex corrections vanish exactly because of the 
mismatch between the angular dependencies of the basis 
function and the order parameters. From group theoretical 
point of view this follows immediately from the observa- 
tion that the gaps belong to B lg , see Eqs. J2]) and (l25j) . 
whereas the basis function transforms according to a dif- 
ferent representation [3D]. With all these we find that the 



This statement is indeed true, and it follows at once that it 
is the F function alone that carries all information about the 
specific physical quantity, its matrix structure in Nambu space, 
the coherence factors involved, and the momentum dependence 
of its kernel. Thus it is obvious that F has to be recomputed 
for each physical quantity. 




w/A 

Fig. 4. (Color online) Zero temperature A\ g Raman 
spectra in dDW+dSC in the unitary limit for a = 0, 
0.025, 0.05, 0.1, 0.2, 0.3, and 0.5, from top to bottom (at 
w = 2A). 



quasiparticle contribution gives the total F function. In 
the continuum limit, where 

A£f 4 (k) = cos(40), (32) 
the following closed analytical formula is obtained 



F(u,u') = 



1 1 



E 



15 u' 2 -u 2 \VT 



= (14- 18u 2 + 20W 



+ 28m 4 - 60uV - 24w b + 40uV) + 



K 



y/i — u- 



:{7u 2 -15uu' 



16u 4 + 40uV + 24u 6 - 40uV) - [« o tt'H , (33) 



where the arguments of the elliptic integrals are the same 
as in Eq. (|19p . Substituting this expressions in the nor- 
mal and pair-breaking integrals in Eq. (|30p yields the fre- 
quency dependence of the susceptibility in the dDW+dSC 
state. At zero temperature it is shown in Fig. 2] for differ- 
ent values of impurity concentration. The spectra do not 
depend on the relative sizes of the dDW and dSC order pa- 
rameters. The only energy scale, apart from the scattering 
rate of course, is the amplitude of the effective d-wave gap 
A. This could be inferred already from Eq. (EHI) . reflect- 
ing the fact that Raman scattering in A\ g symmetry is a 
type I process in both particle-hole and particle-particle 
condensates, see Table [TJ If we take a closer look on the 
plot we can see that there is a very small bump around 
lu = A, which is again caused by transitions from and to 
the impurity band at the Fermi energy, just like in the 
conductivity. What is easy to see though is that with in- 
creasing impurity concentration the peak gets suppressed 
and smoothed out as the lineshape approaches the normal 
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state result [ST] 

lm X t i = 2N ( 1 t\) 2 uj2 ^ 2rr 04) 

Due to the lack of vertex corrections, here the transport 
lifetime equals the total lifetime, which is related to the 
total scattering rate as usual via r -1 = 2F [57] . 

Wrapping up our discussion on the fully symmetric 
representation, for a = the pure result is found to be 

Im X £? 4 = 2iV ( 7 ^ 4 ) 2 tanh (^) 1 

{x- 1 [(7 - 8a; 2 + IQx^Kix) 
-(7-12 a; 2 + 32. T 4 )i?(x)] x<l, 
(ll-28x 2 + 32x 4 )K(x- 1 ) { ' 

-(7-12x 2 + 32x A )E(x- 1 ) x>l, 

with x = u/(2A). In the single component ground states, 
where either A\ = or Ai = 0, it transforms correctly to 
the results obtained in dSC [40] and dDW [5]. 

4.2 Raman spectra in B\ g symmetry 

Next we study Raman scattering in B\ g symmetry. Here, 
the basis function coincides with the angular dependen- 
cies of the dDW and dSC order parameters, and as such 
with that of the full gap as well, and for this reason we 
suspect that this will lead to interference effects between 
the Raman vertex and the gaps [40 54 . Restricting our 
analysis to the the continuum limit as before, where 

Afi' 2 (k) = cos(20), (36) 

the quasiparticle result is again independent of the relative 
sizes of /Ai,2, it reads 

F, p (u, „') = -^— 2 | s/T^ (-nn> + | + £ ) E 

+ 7r^H + 5 + T) A "- [ """'| (37) 

but the vertex corrections do change this situation. The so- 
lution of Eq. (|28p in this symmetry channel is very lengthy 
and we shall only quote here the final result. The technical 
details are given in Appendix [21 The analytically contin- 
ued solution has the form 

ylfi 9 2 (k, u, u') = cos(20)p 3 o- 3 + A(i Pl a 3 ) + B(ia 2 ), (38) 

where the last two terms are the corrections due to dis- 
order. The explicit expressions for the scalar functions 
A(u, u') and B(u, u') are given in Eqs. (g5| and (|g9"|). With 
the impurity renormalized vertex at hand we can now go 
back and evaluate Eq. ([27)1 . After analytically continu- 
ing to real frequencies we find the general expression in 
Eq. (J30J with the total F function 

F(u, u') = F qp (u, u 1 ) + ^AQ + ^-BQ, (39) 




0.5 1 1.5 2 2.5 3 



lu/A 

Fig. 5. (Color online) Zero temperature B\ g Raman spec- 
tra in dDW+dSC in the unitary limit for a = 0, 0.025, 
0.05, 0.1, 0.2, 0.3, and 0.5 from top to bottom (at u = 2 A), 
and A x / A = 2/3. 

where Q{u, u') is given explicitly in Eq. ( |49g[ ) in Appendix. 

This is the point where we can explicitly see what 
we have conjectured in the beginning of this subsection. 
Namely, the interference of the angular dependencies of 
the bare vertex and the order parameters resulted in ad- 
ditional terms in the Raman susceptibility. If the ground 
state is dDW alone, where A 2 — 0, the contribution from 
ladder diagrams is AQ. In dSC the roles simply inter- 
change and we find BQ. In general when both orders con- 
tribute to the coexisting phase neither of these correction 
vanish. The B\ g Raman spectra in the dDW+dSC ground 
state is shown in Fig. [5j We can observe that the peak at 
2A dominates the spectra. This is in fact not that sur- 
prising, because the type I coherence factor is so, as op- 
posed to the conductivity, that it facilitates this transition 
in both the dDW and dSC components of the coexisting 
phase, thus making it more pronounced. With increasing 
impurity concentration the cubic low frequency behavior 
of the pure system changes to linear tendency, the sharp 
peak widens and looses of its intensity as the lineshapes 
approach the normal state result [57 

Im X g^ 2 = 2^(7^) ^ f {2ry (40) 

Apart from trivial substitutions of indices this is precisely 
the same as Eq. (IM)) . showing that in an isotropic metal 
the different symmetry spectra cannot be distinguished. 
There is a remark though we would like to make at this 
point. We see that it is again, interestingly, the total life- 
time of electrons r _1 = 2F showing up in Eq. (j4"0")l . and 
not the transport lifetime, although there are vertex cor- 
rections in this symmetry channel. The reason lies in the 
fact that in our approximation the scattering potential U, 
see Eq. © , does not depend on the momenta of incoming 
and scattered particles. Hence its expansion in the basis 
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functions of the point group consists of only the momen- 
tum independent L = term, which cannot change the 
lifetime in the L = 2 channel [57] . 

Again, we finish this subsection by stating the pure re- 
sult in dDW+dSC, for which closed analytical expression 
exists 



Im X £f 2 = 2iY (7fi 9 2 )^anh 



1 



47V 3 



\ X - 1 ((2 + x 2 )K(x)-2(l + x 2 )E(x)) x<l, 

\(l + 2x 2 )K(x- 1 )~2(l + x 2 )E(x- 1 ) x>l, { ' 

where x = oj/(2A). If we are on the phase diagram where 
the system does not support dDW ordering and only d- 
wave superconductivity exists, we get back correctly the 
known result of dSC [50]. The reverse statement is also 
true: for A — A\ the above expression yields the dDW 
result as well [5]. 



4.3 Raman spectra in B 2g symmetry 

Last but not least let us summarize the results obtained for 
B 2g symmetry. Here the basis function is sin(afca;) sin(a/cj,), 
and the Nambu matrix appearing in the effective density 
is M — (T3, resembling the true charge density. Conse- 
quently, in dDW Raman scattering in this symmetry is 
a type II process, see Table [1] Quite similarly to what 
we have discussed in Ai g symmetry, one finds that the 
vertex is not renormalized by impurities because of the 
mismatch between the basis function and the d-wave or- 
der parameters. They belong to different representations, 
hence the orthogonality. Also, it is worth mentioning that 
though the vertex possesses the same Nambu structure as 
the density itself, still there is no Coulomb screening in- 
volved, because the L = and L = 2 basis functions are 
orthogonal too. This leads to the vanishing of the cross- 
susceptibility ([0-3, 72A2CT3]), so the screening RPA series 
cannot build up. From physical point of view the reason 
is that in case of short wavelength intracell density fluc- 
tuations the Coulomb forces are irrelevant |40) . 

In the continuum limit, where the basis function sim- 
plifies to 

Af!f 2 (k)=sin(20), (42) 

the total F function is again obtained from the one bub- 
ble calculation, from Eq. (|3Tj) . Due to the opposite coher- 
ence factors of dDW and dSC even the quasiparticle result 
shows explicit dependence on the relative sizes of the gaps 



F(u,u') = - 

11' 



1 



u' 2 — u 2 



- u 2 u(u + u')(E - K) 



+ f(^f ) "I --"'^ '-- , "- > ^l ) - I" ~ "'] }■ <«> 



where, again, the argument of K and E is the same as 
in Eq. (fT5]). The spectra in the dDW+dSC ground state 
are shown in Fig. [5J The main qualitative difference here 




w/A 

Fig. 6. (Color online) Zero temperature B 2g Raman spec- 
tra in dDW+dSC in the unitary limit for a = 0, 0.025, 
0.05, 0.1, 0.2, 0.3, and 0.5 from bottom to top (at u = 3A), 
and A x I A = 2/3. 



compared to the other spectra is that the obvious peak 
at 2A is missing. This is again due to the mismatch be- 
tween the angular dependencies of the bare vertex and the 
gap functions, the same effect that lead to the vanishing 
of vertex corrections. With increasing pair-breaking the 
curves smoothly approach the metallic result |57) . just as 
they do in the other representations 



Im X ^ 2 = 2iVo(7r: 9 2 ) 2 



pry 



(44) 



The pure result for a = can again be expressed in 
closed form 



B L 2 J 2 = 2N (^J 2 ) 2 



tanhf^U 
A J \4TJ 3 



x- 1 ((l-x 2 )K(x)-(l-2x 2 )E(x)) x<l, 
\2{l-x 2 )K(x- 1 )-(l-2x 2 )E(x- 1 ) x>l, 1 ' 

with x = lu/(2A). As expected, in a single component 
dSC with A — A% we reobtain the known result and line- 
shape [JO] • However, the other extreme where only dDW 
exists is rather interesting and somewhat controversial as 
far as the literature is concerned. In pure system Eq. (1451) 
reveals that because of the type II behavior there is no 
absorption at any frequency in dDW. This result is com- 
pletely analogue to the optical conductivity of pure su- 
perconductors, which too vanishes in the long wavelength 
limit, again because of the type II coherence factor [13]. To 
further deepen the analogy we find that in the presence of 
non-magnetic impurities, exactly as in the case of the in- 
frared response in superconductors [27141158] , the Raman 
intensity becomes finite, but is certainly very small. This 
result of ours is in precise agreement with the findings of 
Ref. g], which, though focuses on clean cuprate super- 
conductivity, comes to the conclusion that in the strongly 
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underdoped region of the phase diagram the Big Raman 
response is very small due to the relation A2 <C A\. We 
must note at this point, however, that other authors [5| 
find finite Big Raman susceptibility in pure dDW, which 
in our opinion is probably not correct. The discrepancy in 
their calculation originates most probably from the fact 
that they derived the Raman susceptibility from the con- 
ductivity, and might have overlooked the different Nambu 
structure of the current operator (^3) and the effective 
density (po = 1) in this symmetry. 



5 Conclusions 

In summary we have studied coherence effects and dis- 
order in the ground state of a quasi two-dimensional co- 
existing d-wave density wave and ci-wave superconductor 
(dDW+dSC). The main purpose of the paper was to gen- 
eralize and extend available calculations on in-plane op- 
tical conductivity and electronic Raman response in the 
sense to incorporate impurity scattering on a microscopic 
level and to allow for the possibility of a coexisting phase 
with two order parameters, especially relevant in under- 
doped cuprates. 

We calculated the frequency dependence of the in- 
plane optical conductivity and the superfluid density in 
dDW+dSC in the presence of impurities in unitary limit. 
The dDW component of the ground state results in finite 
absorption at all frequencies, even in the pure system, and 
in particular in a strong peak around twice the amplitude 
of the effective d-wave gap, which is reminiscent of the 
type I coherence behavior known from density wave the- 
ory. As to the superfluid density, we found that due to the 
transfer of spectral weight from zero to finite frequency, 
not even in the pure system would it reach unity, instead 
it is reduced according to the the very simple expression: 
I- A\/ A 2 . 

We determined the electronic Raman susceptibility in 
dDW+dSC. The calculations were carried out for the Ai g , 
Big and Bi g symmetries. In accordance with available re- 
sults on Raman lineshapes of cuprates, the spectra are 
strongly channel dependent. The qualitative differences, 
however, continuously disappear with increasing impurity 
concentration as the curves approach the Lorentz-like metal- 
lic form. We found that the dDW and dSC components of 
the coexisting dDW+dSC state contribute practically the 
same in A\ g and B\ g symmetries, because here the coher- 
ence factors in both condensates are type I. As opposed 
to this, in Bi g the Nambu structure of the Raman vertex 
changes from P3U3 to 03, which at the same time changes 
the dDW coherence factor too from type I to type II. This 
in turn leads to the interesting conclusion that the Raman 
intensity in this channel is dominated by the dSC contri- 
bution, the dDW part survives only in dirty materials. 

This work was supported by the Hungarian National Research 
Fund under Grant No. OTKA K 72613. One of us (A. Vanyolos) 
would like to thank the hospitality of the Max-Planck-Institut 
fur Komplexer Systeme in Dresden (MPIPKS), where a part 
of this work was done. 



A Vertex corrections to B\ g Raman vertex 

In this appendix we summarize the technical details of the 
solution of the vertex equation for B\ g symmetry. The cal- 
culation is carried out in the unitary limit. In this sym- 
metry M = P3U3, and the bare vertex in the continuum 
limit takes the form 



AfL 9 2 (k) 



cos(2(/>). 



(46) 



In unitary limit the T-matrix is given by Eq. (|29[) . Direct 
substitution of the bare vertex into the right hand side of 
Eq. (|28[) shows that the /33173 matrix structure cannot be a 
consistent solution alone, because due to the interference 
of cos(20) with the angular dependencies of ^^(k), new 
terms appear, proportional to P1CT3 and 02. Therefore, the 
first iteration already suggests that we have to allow for 
the existence of new terms in Nambu space, and this leads 
us to the ansatz 

4fl 8 2 (k, u, u') = cos(20)p3CT3 + A{i Pl <j 3 ) + B(ia 2 ). (47) 



Substituting this into both sides of Eq. (|28l) one finds that 
this is indeed a consistent solution, and the problem is 
now reduced to finding the solution of two coupled linear 
equations for the unknown coefficients A and B. For these 
in turn we get 



A = 

B = 



be-a(f - 1) 
(c-l)(/-l)-de' 

ad — b(c — 1) 
(c-l)(/-l)-de' 



where 



*r VT^Vi-w 2 Ax , 

a = -tt-t .„„. —Q{u,u ), 



2A 



uu'KK 1 



b = 



rJVl -u 2 \fl^ 



A 

'a. 



2A 



uu'KK' 



Q(u,u'), 



(48a) 
(48b) 

(49a) 
(49b) 



71-r y/i - iivi - w 2 1 fu-u' 
2 A uu'KK' u' - u\u + u' 

( uK u'K' 



VT^u 2 



,/^^HZ£EZ^^Q (uX) , (49 d) 



uu'KK' 



A 2 



e = d, 



(49e) 



nT \/l — u 2 \/ \ — ■ 



1 



2A 



uu'KK' 



u' — u \ u + u' 



uK 



u'K' 



s/T^u 2 



2Q 



A, 



(49f) 
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and 



u + w \ 



y/l - u' 2 



- v 7 ! - u 2 E - 



u 2 K 



(49g) 



In these expressions the argument of K and E, t he ellip- 
tic integrals of the first and second kind, is — u 2 , 

whereas in K' and E' it is l/\/l — u' 2 . 
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